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Abstract. We investigate the dissipativity properties of a class of scalar second order parabolic 
partial differential equations with time-dependent coefficients. We provide explicit condition on the 
drift term which ensure that the relative entropy of one particular orbit with respect to some other 
one decreases to zero. The decay rate is obtained explicitly by the use of a Sobolev logarithmic 
inequality for the associated semigroup, which is derived by an adaptation of Bakry's T— calculus. 
As a byproduct, the systematic method for constructing entropies which we propose here also yields 
the well-known intermediate asymptotics for the heat equation in a very quick way, and without 
having to rescale the original equation. 

Resume. Cet article propose une etude du mecanisme de dissipation d'entropie pour une classe 
d'equations aux derivees partielles paraboliques dont les coefficients dependent du temps. Sous des 
criteres formules explicitement en terme des coefficients, nous etablissons la decroissance exponen- 
tielle de l'entropie relative d'une orbite par rapport a une autre, pour des equations n'admettant pas 
de solution stationnaire. La methode utilisee repose sur l'obtention d'une inegalite de type Sobolev 
logarithmique pour le semi-groupe associe, grace a une adaptation du critere de Bakry-Emery. 
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1 Introduction 

1.1 The Kullback-Leibler Distance as a Particular $-Entropy 

Given two probability densities u, v on M. d , the entropy of u relative to v (also known in information 
theory as their Kullback-Leibler distance (see although it is not a distance) is defined by 

H(u\v) := [ u(x) In ( ^rl) dx. (1) 
Jn* \v{x)J 
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Although this quantity does not satisfy the triangle inequality, it is always nonnegative and vanishes 
only when u = v. These two facts are immediate consequences of the well-known Pinsker inequality 
(see [12]): 

H{u\v) > -\u-v\ 2 Ll . 

Thus the quantity H{u\v) may provide some notion of "distance" between u and v. In Partial 
Differential Equations H(u\v) may be useful in studying the asymptotic behavior of a dissipative 
system. In this context v is in general a stationary solution (or in physical terminology a detailed 
balance equilibrium), u is the orbit of some Kolmogorov (or some other parabolic) equation, and 
H{u\v) is a decreasing function of time. In many cases this fact can be combined with some clever 
inequalities to show that v in fact attracts u in some appropriate metric. This is the basis of the 
well-known entropy dissipation method (see H3|), which has been used to good advantage in many 
examples such as (linear or nonlinear) parabolic equations, kinetic equations, etc... 

As is well-known from information theory and Statistical Physics Q is a particular instance of 
a more general class of entropies: 

H*(u\v):= [ v(x)^(^]dx, (2) 
J Rd \v(x)J 

where $ is any convex function defined on [0, oo[. Formula Q corresponds to the particular choice 
$(z) = zlogz which has an interesting extensivity property jH], but as a general rule the dissipation 
of entropy is a convexity property which has little to do with the specific properties of the z log z 
function. 

In this paper we investigate a class of linear parabolic equations which due to the presence of 
time-dependent coefficients have no stationary solution, but for which orbits still do come together 
in the "metric" given by (£Q) or © . More precisely, we will give explicit conditions on the coefficients 
which ensure that the quantities Q or decrease to zero for large time, with quantitative bounds, 
for any two orbits u and v, that is, even when v is a non-stationary solution. 



1.2 The Entropy Production for Linear Scalar Advection-Diffusion Equations 

Let us consider a general multidimensional linear equation in the form: 



du 

— + divJ(n, Vii) 



0, 



(3) 



where J is some specified flux function. Let $ be any convex function, and assume that J is such 
that this equation preserves positivity. Given two (time-dependent or not) positive solutions u and 
v of this equation, we define the quantity H®(u\v) at any time t by 



H*{u\v) := I v(t, 



x)$ 



u(t, x) 
v(t, x) 



dx 



(From here on, the notation dx will mean d— dimensional Lebesgue measure). The convexity of <1? 
implies a lower bound for H^(u\v ) as follows: using Jensen's inequality with the probability measure 
j V ^ x we obtain the inequality 



v(t, x)& 



u(t, x) 
v(t, x) 



dx > v(t, x) dx& 



( J Rd u(t,x)dx \ 
\f Rd v{t,x)dx)' 



If we assume that $(1) = and that u and v have the same integral (which we only have to assume 
at time zero, since equation ((HJ) will preserve mass for reasonable solutions), we see that H®(u\v) 
remains nonnegative. 

Assuming that u and v vanish at infinity, a straightforward computation then yields: 



£fl*(u|t,)= / WHWH 

dt v 1 ; J Rd \vJ \v 



J(«,Vtt) J(v,Vv) 



u dx. 



(4) 



Perhaps the best known case of this general formula is the case of a one- dimensional Kolmogorov 
(also known as Fokker-Planck equation) equation admitting a detailed balance equilibrium m. In 
this case (JSJ takes the form: 



dt dx\ dx\m. 

Taking v = m in (@J) we obtain : 



^ff*(«|m) = - ! W- 
dt l K d \m 



d / u 
dx \m 



2 

Dm dx, 



in which one recognizes a generalization of the familiar Fisher information, the usual Fisher infor- 
mation corresponding to the case where D(x) = 1 for all x, together with the choice <3?(r) = rlogr. 
Let us now go over to the case of a general linear second-order scalar advection-diffusion equation: 

du 

— + div (b(t,x)u-a(t,x)Vu) = 0. (5) 

Here b is a given vector field, and a is a given diffusion matrix. For any two positive solutions u, v 
to this equation, (@J) now becomes: 



v / V v 



vdx. (6) 



As expected, this formula shows that linear transport does not play any role in the entropy production 
(the expression does not involve the velocity field b): in physical parlance, diffusion here is the only 
irreversible process. Let us emphasize that, from the convexity of and the positive definiteness 
of the diffusion matrix a, we obtain that the relative entropy H(f\g) is a time-decreasing quantity, 
whatever the solutions u, v are, and whatever the coefficients 6, a are. In particular for 
arbitrary time-dependent coefficients b and a the system will have no detailed balance equilibrium , 
i.e. the problem 

b(t,x)u(t,x) = a(t,x)Vu(t,x) 

will have no solution at all. Note that in the context of Markov processes, a similar dissipation 
property was exhibited by Yosida, and Kubo (see [§])■ 

The natural question then arises to investigate under which (sufficient) conditions on the coeffi- 
cients of the equation does the entropy decrease to zero. In the classical setup where one investigates 
the trend toward a stationary solution, it is well-known that for a large class of such stationary so- 
lutions, some Logarithmic Sobolev Inequality is available This fact can be used to obtain a 
Gronwall-type inequality for the entropy, thereby yielding exponential decay. The supplementary 
difficulty here is that the measure relative to which the entropy is computed moves along the flow, 
in such a way that classical conditions which ensure that a logarithmic Sobolev inequality will hold 
cannot be checked a-priori. 

In the next section we revisit a well-known prototype, the Ornstein-Uhlenbeck equation with 
constant drift. This example shows that depending on the nature of the drift the entropy may decay 
to zero in an exponential or algebraic fashion, or converge to a nonzero value. Section 3 collects the 
technical tools needed to show that for the case where the diffusion matrix is the identity matrix, 
the solution of the evolution problem will satisfy the logarithmic Sobolev inequality at all positive 
times. The asymptotic behavior of the entropy is obtained as an easy corollary in Section 4. Finally 
in Section 5 we show that at least for the heat equation (but we believe for a much larger class of 
parabolic problems), the choice of the fundamental solution for v provides a very quick proof of the 
classical Gaussian intermediate scaling. 

2 The fundamental example: the Ornstein-Uhlenbeck process 

Let us consider the simplest case. Denote by (X) t>0 the solution of 



dX t = V2B t - XX t dt, 



where (B) t>0 is a standard Brownian motion on R and A G 1 is a constant. This equation can be 
solved as follows: 

X* = X e- Xt + V2 [ e A ( s -') dB s . 
J o 

As a conclusion, the measure P t (-)(x) which is defined as the law of X t knowing that X = x is the 
Gaussian measure with mean xe~ xt and variance (1 — e~ 2A ')/A. One can then compute the relative 
entropy of P t (-)(y) with respect to P t (-)(x): 

a{t) := H(P t (.)(y)\P t (.)(x)) = ^~^y 

since Pt(-)(y) and Pt(-)(x) have the same variance. Of course, in the case when A = 0, the above 
formula has to be understood as 

H(P t (.)(y)\P t (.)(x)) = 
As a conclusion, three different behaviors can occur: 

• if A > 0, then a decreases exponentially fast to 0, which is natural since Pt(-)(x) converges 
exponentially fast to its invariant measure M(0, 1/A). 

• if A = 0, then a still goes to zero although Pt(-) does not converge to a probability measure, 

• if A < 0, then a converges exponentially fast to a nonzero limit: 



-K*-y? = _K X _ v) 2 _ Kx-v? xt > _K 

2(1 -e xt ) 2 { V) 2(1 -e xt ) t-oo 2 l 



3 The local $-Sobolev inequality for inhomogeneous diffusion Semi- 
groups 

3.1 Notations 

In this section we consider the family of formal elliptic partial differential operators (L t )t>o defined 
by 

d d 

L tf(x) := ^2 a ij {t,x)d ij f{x) + ^2bi{t,x)dif{x), (7) 
i,j=l i=l 

where (ciij(t, -))i<i j<d ^ s a definite positive diffusion matrix and b(t, ■) is a given vector field on 
M. d , defined for all t > 0. Let us suppose that the coefficients are smooth functions of (t,x). This 
family of operators (L t ) t>0 generates a inhomogeneous Markov semigroup which we will denote by 
{Ps,t)o<s<t in the following sense. Writing as usual a as a = era , one can associate to (L t ) t>0 the 
solution of the following SDE: 



A*' r = x + fb(X x s > r ) ds + V2 fa(X x s > r ) dB t 

J r J r 



where (B>t) t>0 is a standard Brownian motion on M. d . Semigroup and probabilistic approaches are 
linked by the fundamental relation 

P 9 , t f(x) :=Ef(Xr). 

The Markov property of X can be translated into a composition rule for the semigroup: for every 
s < t < u, 

P s , u f(x) = E[f(X s u >*)} =E[f(x l u Xr )] = E[P t , u f(Xn) = Ps,tPt,uf(x). 



This semigroup satisfies the well-known Kolmogorov equations: 

d s P s , t f = -L s P s , t f, d t P Sit f = P Sit L t f. (8) 

Let fi be a probability measure on M. d and u the density function of the law of X t knowing that 
C(Xq) = [i. Then, for every smooth function /, 



E(/(X t )) =E(E(/(X t )|X )) = J P 0tt f(x)fi(dx). 
The ltd Formula ensures that, for every smooth function / and s <t, 

f(X t ) - f(X s ) - J L r f(X r ) dr 

is a martingale. In other words, 

Ef(X t ) - Ef(X s ) - y*EL r /(X r ) dr = 0. 

As a consequence, 



u(t,x) — u(s,x) — / L*u(r,x)dr 



f(x) dx = 0, 



and u satisfies Equation © with initial condition /* in a weak sense. 

Following let us associate to L t the two bilinear forms F(t) and defined by: 

m(f,9) ■= \[L t {fg)-gL t f-fL t g\, 
m(f,9) ■= \[L t Y{f,g)-Y{g,L t f)-Y{f,L t g)]. 

We will write T(t)(f) instead of T(t)(J,f) and E (*)(/) instead of T 2 (t)(f, /). 

Remark 3.1 One can c/«ec& that 

1 d 

r(t)(/, 5) = V/ • a(t, -)Vg = -J2 EN + 0^/^- 

i=i i=i 

Remark 3.2 T/ie expression 0/F2 «s much more complicated in the general case. In the simple (but 
informative) case when a(t, •) is the identity matrix, it is very easy to check the following formula: 

r 2 (i)(/) := \\Hess{f)\\l - Vf ■ Jac(b(t))Vf, (9) 

where Hess(-) (resp. Jac(-)) stands for the Hessian (resp. Jacobian) matrix, and \\B\\2 denotes the 
Hilbert- Schmidt norm. 

Remark 3.3 Notice that 

Vf • Jac{b)Vf = Vf ■ SJac(b)Vf, 
where SJac(-) stands for the symmetric part of the Jacobian matrix i.e., 



SJac(b)i 



Jac(b)ij + Jac(b)ji 



The antisymmetric part of the Jacobian of b brings no contribution in our study. One can think 
about the following explicit example: consider the 2- dimensional process X solution of the following 
SDE: 

dX t = dB t -( _^ j J X t dt. 

The antisymmetric part of the drift induces a rotation whereas the symmetric part ensures the con- 
vergence to equilibrium. 



Let $ : I -> M be a smooth convex function denned on a closed interval Z of R not necessarily 
bounded. Let p be a positive measure on a Borel space (O, T). The ^-entropy functional Ent* is 
defined on the set of ^-integrable functions / : (O,^ 7 ) — > (Z,B(Z)) by 

EntJ(/) = jj>(f) dp-<S> (Jj dp)j . 

In what follows, p is a probability measure. As a consequence, J n f dfi £ Z and the definition make 
sense, in the sequel, one has to make an extra assumption in order to derive interesting functional 
inequalities: 

(u,v) i — > <&"{u)v 2 is non negative and convex on Z x Z. (10) 

Remark 3.4 The classical variance and entropy are ^-entropy junctionals respectively associated 
to x i— > x 2 onZ = R and x t— ► x log x on Z = [0, +oo) . 

Definition 3.5 The semigroup (-fs,t)o< s <t * s sa *^ ^° satisfy a local <fr-Sobolev inequality with cons- 
tants (C s ,t) 0<s<t if for all s <t and smooth function f , 

Ent| t (/) := PsAHd)) ~ *(P*,t9) < Cs,tP,,t{*"(f)T(t)(f))' 

Remark 3.6 Under the so-called Bakry-Emery criterion, 

3p e R, V/ smooth, T 2 (/) > pT(f), 

homogeneous diffusion semigroups satisfy a Poincare and a logarithmic Sobolev inequality (see OV). 
As a generalization, Q-Sobolev inequalities can also be established (see 17]). 

Our aim is to take into account the time dependence of the coefficients of the diffusion process. 
We will show that the appropriate adaptation of the Bakry-Emery criterion to that situation is as 
follows: 

3p : p(t) GR, V/ smooth, T 2 (t)(f) + ^d t T(t)(f) > p{t)T(f), (11) 
where dtT(t) is defined as 

d 

d t T(t)(f,g) := dmjitrWdjg. 
».i=i 

The key point in the homogeneous and diffusive case is to get the following commutation relation 
(which turns out to be equivalent to Bakry-Emery criterion): 

In the following subsection we derive such a commutation relation in the inhomogeneous case. 
3.2 The commutation relation 

Let s and t be two fixed times, with < s < t. The key point is the following lemma, which describes 
how the dissipative mechanism tends to flatten gradients: 

Lemma 3.7 Suppose that the family of operators (L t ) t>0 defined in Q) satisfies (|TT]l . For any r 

between and t, the following inequality holds true: 

^T(r)(P T , t g) < exp (- J* p(u) du\ P T , t (V^W^) ■ (12) 
Proof . For all u E [r, t] , we define /3(u) by 

P{u) = P T JyjT(u)(P Uit g)\ 



and compute its derivative by using ©. The crucial assumption that L u is a diffusion operator 
ensures that 

= <S>'(g)L u (g) + <5>" {g)T(u)(g). (13) 

In order to make the exposition clearer we denote P u ,t9 as h and T stands for T(u). A straightforward 
computation leads to: 



13' {u) = P T J L u ^/fP^g + 



2y/F(P u , t g) 



{-2T{P u , t g,L u P Ujt g) + d u T(P U}t g)} 



p 

1 T,U 
P 

1 T,1l 
P 



1 



-L u Th 



1 T(h,Luh) | d u T(h) 



2VTh 4(r7i)§ vTh 2VTh 

2{Th)L u Th - A(Th)T(h, L u h) - TTh + 2(Th)(d u Th) 



4(I7i)5 

4(Th)T 2 (u)(h) - TTh + 2{Th)d u Th 



Therefore we obtain: 

0{u) - p{u)(3{u) = P T , U 



4{Th)i 



4(Th)(T 2 (u)(h) + (l/2)d u Th - p{u)Th) - TTh 
4(I7t)t 



Following one can show, thanks to the diffusion assumption, that the criterion tfTTfl implies that, 
for all smooth functions /, 



r 2 (n)(/) + (i/2)^r/- P (n)r/> 



IT/ 
4T/ ' 



Remark 3.8 In the case when a is the identity matrix, one can easily get derive this inequality from 
the criterion and the Cauchy-Schwarz inequality. 

As a conclusion, f3 satisfy the following differential inequality: 

0{u) > p(u)P(u). 

In other words, the function 

u i — > (3(u) exp ^— J p(v) dv 
is an increasing function on the interval [r, t] which implies that 



/5(r) exp 
This is precisely the desired inequality. 



p(u) du I . 



□ 



3.3 Local $-Sobolev inequalities 

Theorem 3.9 Suppose that the family of operators (L t ) t>0 defined in satisfies ifTTjl . Then for 
any times s,t with < s < t and any positive function g, P s j satisfies the following Q-Sobolev 
inequality: 

Ent% t (g)<c(s,t)P s>t ($"(g)T(t)(g)), 
where the constant c(s, t) can be chosen as: 

c(s,t) = / exp I —2 / p(u)du I dr. 



Proof. Consider the function a : [s, t] — ► K denned by : 

a(r) := P s , T (*(P r ^)), 

Let us compute the derivative of a: 

a'(r) = P SjT {L T {${P Tjt g)) - &(P Ttt g)L T P Ttt g) 
Thanks to the fact the diffusion assumption, (|T3|l ensures that 

a'{r)=P s ^"{P T>t g)T{T){P T>t g)). 
The commutation relation lfT2|l ensures that 

<f>"(P Tjt g)T(T)(P Tjt g) = ^\P Ttt g){jT{T){P T ^ 2 

< exp ^-2 J p(u) dv) $"(P Ttt g)P T)t (yr(t)(g) 

Jensen inequality with the bivariate function (u,v) \— > <&"(u)v 2 (which is assumed to be convex 
according to ifTHjl ) ensures that 



$\p T , tg )p T , t ^/fWg)) < Pr,t{^"{gW){g))- 

As a conclusion, 

a'(r) < exp f-2 J p(u) dv\ P s>T P T , t (^"(g)T(t)(g)) = exp (-2 j\(u) dv\p s>t ($" \g)T(t)(g)) . 

Since a(t) = P s> t($>(g)) and a(s) = $>(P St t(g)), the result follows upon integration of this inequality 
between s and t. □ 

In the special case of $ : x t— > xlogx, the local logarithmic Sobolev inequality for the semigroup 
(P Si t) 0<s< £ can be stated as follows: 

Corollary 3.10 Suppose that the family of operators (Lt) t>0 defined in satisfies ifTTIl . Then for 
any times s, t with < s < t and any positive function g, P Sj t satisfies the following logarithmic 
Sobolev inequality: 

Ent ftjW == P., l(9 lo g9 , - (P M P„ tg) < o(M)P„ t (^), (14) 
where the constant c(s, t) can be chosen as: 

c(s,t) = J exp ^—2 J p{u)du^j dr. 

Remark 3.11 If for every x G W 1 , the matrix (aij(t,x))ij is bounded by the identity matrix (as 
symmetric bilinear forms), then T(t)(g) < |Vg| 2 and P s> t satisfies the classical logarithmic Sobolev 
inequality: 



Ent Ps Jg) < c(s,t)P t 



|V 5 | 2 



g j 

Remark 3.12 In the case where p(t) = p for all t > 0, we recover the constant 

cM) = f (l-e- 2 ^)) 
provided by J^j in the homogeneous case. 



4 The Kullback-Leibler distance of two orbits of a parabolic prob- 
lem 



We now consider the following parabolic equation © 

dv 

at 



()V (t,x) = L*v(t,x), t>0,xeR d , (15) 



v(0,x) = vq(x), xeR d , (16) 

where L\ stands for the adjoint of L t with respect to the Lebesgue measure on M. d . 

If we assume that the initial data vq satisfies a logarithmic Sobolev inequality, the result of the 
previous section may be used to show that the inequality is propagated in time: 

Theorem 4.1 Assume the initial data Vo satisfies the following logarithmic Sobolev inequality: for 
all smooth functions f , 



f log fv - / fv log / fv < d { 



./r(o)(/H. 



// the family (Lt) t>0 satisfies the criterion (|11|1. Then for any positive time, the solution v(t, •) of 
f75))J7^)) satisfies the following logarithmic Sobolev inequality 



f log fv(t,-)- fv(t,-)log fv(t,-)<d(t) /r(t)(/)«(t,-), 



where 



d{t) := d exp ^-2 J p{r) dr^j + J exp ^-2 J p(r) dr^j dr. (17) 
Proof . For any positive function g we have: 

•,(,) .o g9W „(M) *, = /e[/(x-)] „ oW = / P„,,( 9 lo g9 )(,)„(0„) 



The integrand may be bounded from above by using the local inequality (|14|) with s = 0, and the 
logarithmic Sobolev inequality for vq: 

J P ,t{9 log g)v(x,0)dx < J P 0>t g log P 0>t gv o (x)dx + c(0,t) / Po,t( ^y^ Jv (x) dx 

< d f^^v (x)dx + c(0,t) J P , t (^y (x)dx. 



The first integral in the last inequality may be estimated by 1|T2|) . and this completes the proof. □ 

We are now in a position to estimate the Kullback-Leibler distance between two arbitrary orbits 
of (IS): 

Theorem 4.2 Under the same assumptions as in previous theorem, let u be another solution of 
4J.5)) (i.e., corresponding to different initial data uq). Assume that uq and Vq are positive. Define 
the relative entropy of u with respect to v at any positive time t by 

H(u(t)\v(t)) := / u (t,x)\og^\dx. 

J v(t,x) 

This quantity is then bounded as follows for all positive times: 

H(u(t)\v(t)) < H(u(0)\v(0))c(t), 
where t 

cW = exp ("i W) ds 

and d(t) is the constant defined in i 17\) . 



Proof . The proof is a straightforward application of Gronwall's lemma. Let us set g := ^; from JgJ) 
with <3?(,z) = zlogz, we obtain 



|i/ Wt ) W t)) = -/^»dx. 



Therefore, Theorem 14. II gives the following control: 

jH{u(t)\v{t)) < --L H (u(t)\v(t)), 

which gives the result. Let us conclude this section by indicating what the obtained rate is when 
the quantity p{t) may be taken to be a fixed constant p. If p = 0, then (|T7]l gives d(t) = do + t, 
therefore we get the algebraic decay 

For p / 0, we have 

d(t) = d e^ t + ^(l-e-^ t ), 
2p 

and the integral in lfT%|) may be computed to yield 

2pd e- 2 P t 



c(t) 



1 + (2pd - l)e- 2 ^ ' 

□ 



5 An application to intermediate asymptotics 

In the case where a Kolmogorov equation has a rather trivial (e.g. constant) asymptotic state, it is 
often the case that (due to the self similarity of the underlying Markov process) some appropriate 
rescaling of the orbit shows structure (e.g. Gaussian), a phenomenon termed intermediate asymp- 
totics by Barenblatt pj. To illustrate this point let us consider the one dimensional linear heat 
equation on the entire line: 

du 1 d 2 u m 
— = -— . iel. 
dt 2 dt 2 ' 

It is very easy to check that for large times u converges point wise to zero, but it is also very easy 
to read off from the explicit form of the solution (as given in terms of the heat kernel) the following 
point wise convergence: 

Vtu(yVt,t)^Cexp(-V-), (19) 



where the constant C may be determined from mass conservation. As is well-known, this convergence 
may be obtained by rescaling the equation and constructing an entropy functional for the rescaled 
equation. Let us briefly recall this argument. The first step consists in rescaling the function u by 
setting: 



u(t,x) = a(t)v\—^-,T(t) 

where the scaling functions a(t), (3(t), r(t) are to be chosen so as to make the equation for v as simple 
as possible, while preserving the mass constraint: 



u(t, x) dx = const. 
These two requirements lead to the choice 



v(y, kit) = Viu(yVi,t), 



(20) 



thus to the equation 

dv Id/ dv\ 
dr 2 dy \ dy J ' 

The detailed balance equilibria for this equation are exactly all multiples of the standard Gaussian 
density, i.e. take the form: 

Voo (y) = Cexp f-y), (21) 
for some constant C. By the result of previous section, the relative entropy 

H(v\ Voo ) = / vln — dy (22) 

decreases to zero for large times, which implies that v converges to a fixed point of the form l(2~T1l for 
some C. 

This result, when rephrased in terms of u, is exactly the intermediate asymptotics lfT?H), Note 
that in this case the constant C is uniquely determined from the mass conservation relation: 

u{t,x)dx = / u(x,0)dx. 
Jr 

The key point is now that by using the change of variable x = yy/i in (|22|) one obtains: 

f uit X) 

H(v\v 00 )= / u(t,x)]n— — , dx, 

Jr C^=exp(—^) 

which is exactly the relative entropy of u with respect to the fundamental solution of the heat equa- 
tion. The fact that this entropy is dissipated combined with Pinsker's inequality now immediately 
leads to the intermediate asymptotic without having to resort to any rescaling of the function 
u. 

In other words, the fundamental solution encodes the intermediate asymptotics. 
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